Simulation of Quantum Computation: A deterministic event-based approach* 
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We demonstrate that locally connected networks of machines that have primitive learning capa- 
bilities can be used to perform a deterministic, event-based simulation of quantum computation. We 
present simulation results for basic quantum operations such as the Hadamard and the controlled- 
NOT gate, and for seven-qubit quantum networks that implement Shor's numbering factoring algo- 
rithm. 
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I. INTRODUCTION 

Recent advances in nanotechnology are paving the way to attain control over individual microscopic objects 1-5 . The 
' ability to prepare, manipulate, couple and measure single quantum systems is essential for quantum computation 6 . 
, Candidate systems to implement a quantum computer are ions or atoms in electromagnetic or optical traps, photons in 
cavities, nuclear spins on molecules, quantum dots and superconductors. These technological developments facilitate 
the study of single quantum systems at the level of individual events. Such experiments address the most fundamental 
aspects of quantum theory. Indeed, quantum theory gives us only a recipe to compute the frequencies for observing 
events. It does not describe individual events, such as the arrival of a single electron at a particular position on 
the detection screen 7-10 . Reconciling the mathematical formalism (that does not describe single events) with the 
' experimental fact that each observation yields a definite outcome is often referred to as the quantum measurement 
, paradox. This is the fundamental problem in the foundation of quantum theory 7,8,11 . 

In view of this, it is not a surprise that some of the most fundamental experiments in quantum physics have not 
been simulated in the event-by-event manner in which the experimental observations are actually recorded 12 . One of 
■ the examples are the two-slit experiments in which the individual electron counts build up the interference pattern 9 . 
fS I Other examples are the single-photon beam splitter and Mach-Zehnder interferometer experiments 13 . Within the 
, 1 1 , standard formalism of quantum theory, no algorithm has been found to perform an event-based simulation of definite 
£^ ' individual outcomes in quantum experiments 8 . 

In this paper, we take the point of view that a physical theory such as quantum theory, is a specification of an 
algorithm to compute numbers that can be compared to experimental data 14 . Thinking in terms of algorithms opens 
new possibilities to simulate phenomena for which a proper physical theory is not (yet) available. As discussed 
before, quantum theory is unable to describe individual events in an experiment but it provides an algorithm to 
compute the final, collective outcome. We have already demonstrated that it is possible to construct deterministic 
processes that generate events at a rate that agrees with the quantum mechanical probability distribution, without 
using quantum theory 15 . In this paper we apply these concepts to quantum computation. We present results of 
event-based simulations of single-qubit quantum interference (Mach-Zchnder interferometer), a two-qubit quantum 
circuit (controlled-NOT gate), and Shor's algorithm 16 to factorize N — 15 on a seven qubit quantum computer. 

The event-based simulation method that we describe in this paper is not a proposal for another interpretation of 
quantum mechanics. To avoid misunderstandings, we emphasize that our approach is not an extension of quantum 
theory. We simulate quantum systems without making use of the rules (algorithms) quantum theory provides us. 
However, the final, collective results of our event-by-event deterministic, causal learning processes are in perfect 
agreement with the probability distributions of quantum theory 11 . The event-based simulations build up the final 
outcome event-by-event, just like in real experiments. In physics terminology the entire approach is particle-like and 
satisfies Einstein's criteria of realism and causality 8 . In this sense, our approach can be viewed as a recipe to construct 
event-based systems, that is "classical" models, that behave as if they were "quantum mechanical" . 

The paper is organized as follows. In Section II we briefly recall some basic elements of quantum computation. We 
describe the Hadamard, controlled-NOT (CNOT) and Toffoli gate, involving one, two or three qubits, repectively. 
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We discuss the quantum network for the Mach-Zchnder interferometer and for Shor's quantum algorithm to factorize 
N = 15 on a seven-qubit quantum computer. 

In Section III we introduce the deterministic learning machine (DLM) that is at the core of the event-based 
simulation approach. We present a mathematical analysis that proves that these machines generate events at a rate 
that agrees with the corresponding quantum mechanical probabilities. One of the essential features of (networks of) 
DLMs is that they process one event at a time. After applying a deterministic decision process to the input event 
and sending out an output event, a new input event can be processed. Another essential feature of DLMs is that they 
do not store information about individual events. Networks of DLMs are capable of unsupervised learning 15 but they 
have very little in common with neural networks 17 . The sequence of events that is generated by a DLM is stricly 
deterministic. This is modified in the stochastic learning machine (SLM). The event-by-event learning processes in a 
SLM are still deterministic and causal but the output events are randomly distributed. This modification is necessary 
if we want to mimic the apparent random order in which quantum events are detected in experiments 9 ' 13 . 

In Section IV we first describe the construction of DLM networks and present results from the event-based simulation 
of the Hadamard gate and the Mach-Zchnder interferometer, the latter showing that DLM-based networks correctly 
reproduce quantum interference phenomena. Then we describe how to simulate a CNOT gate using DLM and SLM 
networks. Finally, we present the simulation results of Shor's number factoring algorithm 16 implemented on DLM 
and SLM networks. A summary and outlook is given in Section V. 



II. QUANTUM COMPUTATION 

This section summarizes those aspects of quantum computation that are necessary to understand the examples 
of quantum algorithms that we use in Section IV to demonstrate that local, causal and deterministic processes can 
simulate quantum computers on an event-by-event basis. 



A. Preliminaries 



The state of an elementary storage unit of a quantum computer, the quantum bit or qubit, is described by a 
two-dimensional vector of Euclidean length one. Denoting two orthogonal basis vectors of the two-dimensional vector 
space by |0) and |1), the state |$) of the qubit can be written as a linear superposition of the basis states |0) and |1): 

|$)=ao|0)+ai|l), (1) 

where cto and a\ are complex numbers such that | do | 2 + 1 ai | 2 = 1. The appearance of complex numbers suggests that one 
qubit can contain an infinite amount of information. However, it is impossible to retrieve all this information 10 ' 18,19 . 
The result of inquiring about the state of the qubit, that is the outcome of a measurement, is either or 1. The 
frequency of obtaining (1) can be estimated by repeated measurement of the same state of the qubits and is given 
by |a | 2 (| fll | 2 ) 10 ' 18 ' 19 . 

According to quantum theory 20 , the internal state of a quantum computer with L qubits is described by a unit 
vector (state vector) in a D = 2 L dimensional space (of complex numbers) 6 . Adopting the convention of quantum 
computation literature 6 , the state of an L-qubit quantum computer is represented by 

|$) = a(0...00)|0...00)+o(0...01)|0...01) + ... + o(1...10)|1...10)+a(l...ll)|l...ll) 

= a |0) + ai|l) + . . . + a 2 L_ 2 \2 L - 2) + a 2 i,_ 1 |2 L - 1), (2) 

where in the last line of Eq. (2), the binary representation of the integers 0, ...,2 L_1 was used to denote |0) = 
|0 . . . 00), . . . , \2 L — 1) = 1 1 ... 11) and ao = a(0 . . . 00), . . . , a 2 L_ 1 = a(l . . . 11). As usual, we normalize the state 
vector, that is ($|$) = 1, by rescaling the complex- valued amplitudes a, according to 
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(Hf = L ( 3 ) 



The internal state of the quantum computer evolves in time according to a sequence of unitary transformations 6 . 
A quantum algorithm is a sequence of such unitary operations. Of course, not every sequence corresponds to a 
meaningful computation. Furthermore, if a quantum algorithm cannot exploit the fact that the intermediate state of 
the quantum computer is described by a linear superposition of basis vectors, it will not be faster than its classical 
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FIG. 1: Graphical representation of some of the basic gates used in quantum computation; (a) Hadamard gate; (b) Rotation by 
tt/2 about the x-axis; (c) Rotation by tt/2 about the j/-axis; (d) Single qubit phase shift by <fi; (e) CNOT gate; (f) Controlled 
phase shift by <j> = ir/k; (g) Toffoli gate. The horizontal lines denote the qubits involved in the quantum operations. The dots 
and crosses denote the control and target qubits, respectively. 



counterpart. As the unitary transformation may change all amplitudes simultaneously, a quantum computer is a 
massively parallel machine 6 , at least in theory. 

It has been shown that an arbitrary unitary operation can be written as a sequence of single qubit operations 
and the controlled-NOT (CNOT) operation on two qubits 6,21 . Therefore, single-qubit operations and the CNOT 
operation are sufficient to construct a universal quantum computer 6 . The final state of the quantum computer can 
be calculated by performing the (unitary) matrix-vector multiplications that correspond to the application of this 
sequence (determined by the quantum algorithm) of single-qubit and CNOT operations. According to quantum 
theory, the squares of the absolute values of the elements of the state vector are the probabilities for observing the 
quantum computer in one of its 2 L states. 

The unitary time evolution of the internal state of the quantum computer is interrupted at the point where we inquire 
about the value of the qubits, that is as soon as we perform a measurement on the qubits. If we perform the readout 
operation on a qubit, we get a definite answer, either or 1, and the information encoded in the superposition is lost. 
The process of measurement cannot be described by a unitary transformation 8 ' 10 . Therefore, we do not consider it 
to be part of a quantum algorithm. 



B. Single-qubit operations 

In general, a qubit can be represented by a spin- 1/2 system. The state |$) of a qubit (see Eq.(l)) can therefore 
also be written as a linear combination of the spin-up and spin-down states 10 ' 18,19 : 

|*) = ao|T> + oiU), (4) 

where 6 

|o) = IT) = (J) , |i> = U> = (J)- (5) 

The three components of the spin-1/2 operator S = (S x , S v , S z ) are defined (in units such that K = 1) by 10,18 ' 19 

*-\{ *--\{ ;.;). 

and are chosen such that | f) and | J.) are eigenstates of S z with eigenvalues +1/2 and -1/2, respectively. 
The expectation values of the three components of the qubits are defined as 



(Q a ) = 1/2 - (S a ) , a = x,y,z, 



(7) 
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FIG. 2: Left: Diagram of the Mach-Zehnder interferometer, consisting of two beam splitters, two mirrors, and a device that 
changes the phase of the light wave by <j>. Right: Diagram of the equivalent quantum network. 



where (A) = A qubit is in the state |0) or |1) if (Q z ) = or (Q z ) = 1, respectively. 

A rotation of the state Eq.(4) about a vector v corresponds to the unitary matrix 
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where 1 denotes the unit matrix and v = ^Jv% + v 2 + v 2 z is the length of the vector v. 

For later reference, it is useful to list a few special cases of Eq.(8). The Hadamard operation H and rotations X 
and Y of the state vector by ir/2 about the x and y-axis, respectively, are defined by 6 
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We also introduce the symbol 
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to represent the single-qubit phase-shift operation by a phase <f>. The graphical symbols of H, X, Y, and R(<j>) are 
shown in Fig. 1. The inverse of a unitary operation U is denoted by U. 

The Mach-Zchndcr interferometer is a simple but nontrivial example of a single-qubit system in which the infor- 
mation contained in the phase of the wave function is essential 6 ' 13,22 ' 23 . The schematic layout of the apparatus is 
shown in Fig. 2. N photons enter the first beam splitter through the input channels or 1. The beam splitter 
distributes the photons over its two output channels (0 or 1): The number of photons in these channels is A^o and 
Ni , respectively. If there are no photons in one of the input channels, the beam splitter equally divides the photons 
over its output channels, that is No — Ni — N/2. The photons then propagate and experience a phase shift of 
if they left the beam splitter via channel 1. The photons are collected at another beam splitter. Finally, detectors 
count the number of photons in the two output channels and 1 of the second beam splitter. The number of photons 
in these channels is denoted by N2 and N3, respectively. We assume that no photons are lost in this process so that 

^ = ^0 + ^1=^2 + ^3. 

The relation between input and output amplitudes of the Mach-Zehnder interferometer is given by 6 ' 13,22 ' 23 
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where (ao, a\) and (&o> bi) are the complex-valued amplitudes at the input and output, respectively. The X operations 
represent the beam splitters while the phase shift R((f>) mimics the effect of changing the optical path length. 



FIG. 3: Quantum circuit representation of a CNOT gate where the second (bottom) qubit acts as control qubit and the first 
(top) qubit is the target qubit. 



Let us consider the case where all the photons enter the interferometer through channel 0. Assuming that the 
photons originate from a coherent source, we have (ao,&i) = (cos-0o + i sin -00,0). From Eq.(ll) it follows that the 
probabilities for observing a photon at one of the output channels is given by 

|6 | 2 = sin 2 ^ N 2 = cos 2 ^. (12) 



C. Two qubits: CNOT operation and controlled phase shift 

Computation requires some form of communication between the qubits. It has been shown that any form of 
communication between qubits can be reduced to a combination of single-qubit operations and the CNOT operation 
on two qubits 6,21 . By defmtion, the CNOT gate flips the target qubit if the control qubit is in the state |1) 6 . If we 
take the first qubit (that is the least significant bit in the binary notation of an integer) as the control qubit, we have 

CNOT 21 |$) = CNOT 21 (a |00) +ai|01) +a 2 |10) +a 3 |ll)) 

= ao|00)+o 3 |01) + a 2 |10)+oi|ll) 

= a |0) +a 3 |l) +a 2 |2) +oi|3) 

= ao|0)i|0)2 + a 3 |l)i|0)2 + 02|0)i|l)2 + oi|l)i|l)2. (13) 

where a , . . . , a 3 are the probability amplitudes of the four different states and |0)j and |l)j represent the |0) and |1) 
state of the i-th qubit, respectively. The graphical symbol of the CNOT operation is shown in Fig. le. The dot (cross) 
denotes the control (target) qubit. A CNOT gate in which the control and target qubit are interchanged can be built 
from four Hadamard gates and one CNOT gate 6 . The quantum circuit for this "reversed" CNOT gate is shown in 
Fig. 3. 

The CNOT operation is a special case of the controlled phase shift operation Rji (</>). The controlled phase shift 
operation with control qubit 1 and target qubit 2 reads 



R2i(<p)m = J22i(^)(a |00) + oi|01) + O2|10) + o 3 |ll)) ) 

= aolOO) + d + ^i + (l-^3 |01) + fl2 | 10) + d-e^ + d + e^a,,^ (M) 

Graphically, the controlled phase shift Rji{4> — n/k) is represented by a vertical line connecting a dot (control bit) 
and a box denoting a single qubit phase shift by n/k (see Fig. If). 
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FIG. 4: Quantum network for Shor's quantum algorithm to find the factors of the number N = 15 for the case a — 11 (black 
colored CNOT gates only) and a — 7 (gray colored CNOT and Toffoli gates only) 5 . H denotes the Walsh-Hadamard transform. 
The operations in the dashed box perform a 3-qubit Fourier transform 6 . The operations "2" and "4" perform controlled phase 
shifts with angles n/2 and it/4, respectively. The other gates perform two-qubit (CNOT) or three-qubit (Toffoli) operations. 
The initial state of the seven qubits of the quantum computer is (0,0,0,0,0,0,1). 



D. Three qubits: Toffoli gate 



The Toffoli gate is a generalization of the CNOT gate in the sense that it has two control qubits and one target 
qubit 6 ' 24 . The target qubit flips if and only if the two control qubits are set. Symbolically the Toffoli gate is represented 
by a vertical line connecting two dots (control bits) and one cross (target bit), as shown in Fig. lg. 



E. Seven qubits: Factoring N — 15 using Shor's algorithm 

We now consider the problem of factoring integers. For the case N — 15, an experimental realization of this 
quantum algorithm on a scvcn-qubit NMR quantum computer is described in Refs. 25 and 26. The theory behind 
Shor's algorithm has been discussed at great length elsewhere 6,16,27 . Therefore, we only recall the basic elements of 
Shor's algorithm and focus on the implementation of the algorithm for the case N — 15. The theory in this section 
closely follows Ref. 25. 

Shor's algorithm is based on the fact that the factors p and q of an integer N = pq can be deduced from the period 
M of the function f(j) = a^modN for j = 0, . . . , 2™ — 1 where N < 2". Here a < N is a random number that has 
no common factors with N. Once M has been determined, at least one factor of N can be found by computing the 
greatest common divisor (g.c.d.) of N and a M/>2 ± 1. 

For N — 15, the calculation of the modular exponentiation a^modN is almost trivial. Using the binary represen- 
tation of j we can write a j modN = a 2 " -1 - 7 ™- 1 . . . a 2 - 7 ' 1 a J '° modiV = (a 2n ~ 1 - 5n - 1 mod./V) . . . (a 2jl modiV)(a J °mod./V)modiV, 
showing that we only need to implement (a 2 • J '=modiV). For N = 15 the allowed values for a are a — 2,4, 7, 8, 11, 13, 14. 
If we pick a = 2,7,8,13 then a^modiV = 1 for all k > 1. For the remaining cases we have a 2 *modA f = 1 for all 
k > 0. Thus, for N = 15 only two (not four) qubits are sufficient to obtain the period of f{j) = ajmodN 25 . As a 
matter of fact, this analysis provides enough information to deduce the factors of N = 15 using Shor's procedure so 
that no further computation is necessary. Nontrivial quantum operations are required if we decide to use three (or 
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more) qubits to determine the period of f(j) — modN 25 . Following Ref. 25, we will consider a seven-qubit quantum 
computer with four qubits to hold f(J) and three qubits to perform the Fourier transform to determine the period 
M. 

The quantum circuit that implements Shor's algorithm that factors N = 15 using a = 7 and a = 11 is depicted in 
Fig. 4 25 . Qubits 1 to 3 and 4 to 7 are used as registers to represent j and f(j) = a^modN, respectively. Here, qubits 
3 and 7 are the least significant qubits of these two registers. The initial state of the qubits is (0, 0, 0, 0, 0, 0, 1), that 
is, all qubits are prepared in the state |0) except for qubit 7 which is prepared in the state |1). 

The quantum networks to compute a 3 modl5 for j = 0, . . . , 7 and a fixed input a are easy to construct. Examples 
for a — 7 (six CNOT and two Toffoli gates) and a = 11 (two CNOT gates) are included in Fig. 4 25 . For example, 
consider the case a = 11 = 1 1011) . If j is odd then lPmodl5 = 11 and the network should leave 1 101 1) unchanged. 
Otherwise, lPmodl5 = 1 and hence it should return |0001). The network for this operation consists of two CNOT 
gates that have as control qubit, the same least-significant qubit (qubit 3 in Fig. 4) of the three qubits that are input 
to the Fourier transform. The sequence of CNOT and Toffoli gates that performs similar operations for the other 
cases can be found in the same manner. In the NMR implementation, additional simplications of the a = 7 circuit 
where necessary to render the experiment feasible 25 . There is no need to do this here, so we use the circuits as shown 
in Fig. 4. Elsewhere, we describe a quantum computer emulator (QCE) that simulates models of ideal and physically 
realizable quantum computers 28-30 . The software distribution of QCE 31 contains an implcmcntion of the circuits 
shown in Fig. 4, demonstrating that these circuits work correctly. 

For the quantum network of Fig. 4, we can determine the period M of the function f(j) = a^modN from the 
expectation values of the first three qubits. The state of the quantum computer before it starts performing the 
Fourier transform can be written as 

1 JV-l ( M-1 2M-1 

= 5>'>l/(i)> = 7^1 E \i)\fU))+ E + - 

j=0 ViV [ j=0 j=M 

M-1 

= 7^E(iJ'> + tf +M > + -")i/tf)>. ( 15 ) 

V 1=0 

where, in the last step, we used the periodicity of /(j). Using the Fourier representation of \j) we obtain 

N-l JV-1M-1 
-j=^\j)\f(j)) = ^E E^^ '(l + e^/* ' + e ^M, N ■ + _^ + e 2,ikM { L-l )/ N^ mm 
V 3=0 fc=0 j=0 

1 JV-1L-1 

+ E E e 2 **SI N e™ kML ' N \k)\m), (16) 

where L = [N/M\ denotes the largest integer L such that ML < N. In simple terms, L is the number of times the 
period M fits into the interval [0, N — 1]. The probability p q (M) to observe the quantum computer in the state \q) is 
given by the expectation value of the (projection) operator Q = \q)(q\. With the restriction on f(j) that f(j) = f(f) 
implies j = j' , we find 

m\- (M\ - M f M*qML/N) \ 2 N-MLsm(ngM(2L + l)/N) 

W-p q (M) - N 2 [ sin{nqM/N) ) + N 2 S m(irqM/N) ' ( <} 

The results for p q (M) in the case N = 8 (three qubits) are given in Table I. From Table I it follows directly that 
the expectation values of the qubits are ((Qf) = {QD — (Qf) = 0) if the period M = 1, ({Q\) = (Qf) = 0, 
(Q§) = 0.5) if the period M = 2, ((Qf) = 0.5, (Qf) = 0.375, = 0.34375) if the period M = 3, and ((Qf) = 0, 
(Q|) = (Q|) = 0.5) if the period M = 4. 

Thus, in this simple case of N — 15, the periodicity of f(j) can be unambiguously determined from the expectation 
values of the individual qubits. For a = 7 we find ((Qf ) = 0, (Qf) = 0.5, (Qf) = 0.5) and hence the period M = 4, 
yielding the correct factors g.c.d.(7 2 ± 1,15) = 3, 5 of N = 15. Similarly, for a = 11 we find ((Qf) = 0, (Qf) = 0, 
(Qf) = 0.5) corresponding to the period M — 2 and the factors g.c.d.(ll ± 1, 15) = 3, 5. 



III. EVENT-BASED SIMULATION OF QUANTUM PHENOMENA 



The conventional approach for simulating quantum systems (and quantum computers in particular) is to solve 
the time-dependent Schrodingcr equation of the corresponding system 32 . In practice, this amounts to multiplying 
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TABLE I: Probability p q (M) to observe the state \q) after performing the quantum Fourier transform on the periodic function 
fU) = fti + M) for j = 0, . . . , 7. 



q Pg(M = l) P g (M = 2) P g (M = 3) Pq (M = 4) 

1 0.5 0.34375 0.25 

1 0.0 0.01451 0.00 

2 0.0 0.06250 0.25 

3 0.0 0.23549 0.00 

4 0.5 0.31250 0.25 

5 0.0 0.23549 0.00 

6 0.0 0.06250 0.25 

7 0.0 0.01451 0.00 



the state vector by a sequence of (many) unitary matrices 28-30 . According to quantum theory, the squares of the 
amplitudes of the final state of the quantum computer yield the probabilities for observing the quantum computer in 
a particular basis state. Once these probabilities are known, it is trivial to construct a random process that generates 
events according to these probabilities. 

The approach we propose in this paper is radically different. We construct processes that generate events of 
which the ratio of occurrence agrees with quantum theory. Thus, adopting this method, we don't use concepts of 
quantum theory at all: we don't use wave functions or the Schrodinger equation and do not run into the fundamental 
measurement paradox 8 . In our approach, quantum mechanical behavior is the result of a causal, deterministic (or 
stochastic), event-based process. 

In quantum physics an event corresponds to the detection of a photon, electron, etc. In our simulation approach 
an event is very much the same thing: It is the arrival of a message at the input port of a processing unit. In this 
paper we only consider networks of processing units in which only one message is traveling through the network at 
any time. Thus, the network receives an event at one of its inputs, processes the event and delivers the processed 
message through one of its output channels. After delivering this message the network can accept a new input event. 

The key feature of our approach is a processing unit that we call deterministic learning machine (DLM). A DLM 
is a machine with an internal state that is updated according to a very simple, deterministic algorithm. A DLM 
responds to the input event by choosing from all possible alternatives, the internal state that minimizes the error 
between the input and the internal state itself. This deterministic decision process determines which type of event 
will be generated as output by the DLM. Furthermore, the same process generates different events in such a way that 
the number of each type of event is proportional to the corresponding probabilities of the quantum mechanical device 
that we want to simulate. The message contains information about the decision the DLM took while updating its 
internal state and, depending on the application, also contains other data that the DLM can provide. By updating its 
internal state, the DLM "learns" about the input events it receives and by generating new events carrying messages, 
it tells its environment about what it has learned. This primitive learning capability is the essence of our approach. 

A. Deterministic learning machines 

An extensive discussion of the internal operation and dynamic behavior of a DLM can be found in Refs. 15,33 . 
Therefore, in this paper, we briefly recall the basic ideas. As an example, we take a DLM that we use to simulate 
the Hadamard operation (see Section II) . We first consider a DLM that accepts as input two different types of events 
(0 or 1), each event carrying a message consisting of two real numbers. The DLM learns from the input event by 
updating its internal vector. The internal state of the DLM is represented by a unit vector of four real numbers 
x = (xq, xi,X2,Xs). The first (last) two elements of x are used to learn about the message carried by (1) events. As 
an input event is either of type or 1, the DLM receives a message with two real numbers, not with four. Therefore, 
the DLM uses its internal vector to supply the two missing real numbers. Thus, for a event carrying the message 
yo = (yo,yi), the input vector is v = (2/0, yi, X2, £3) and for a 1 event carrying the message yi = (2/2,2/3)5 the input 
vector is v = (xo,Xi, 2/2,2/3)- Consider now the second possibility in which the DLM receives input in the form of four 
real numbers v = (vo, v\, V2, V3). Then, it is clear that it can skip the step of supplying the missing information and 
use the input v directly. 
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The DLM updates its internal vector by selecting from the eight candidate update rules {j = 0, 1, 2, 3; Sj = ±1} 

w it j = s j\J 1 + a 2 (x 2 - l)Sij + axi(l - S itj ) (18) 
the update rule that minimizes the cost 

C = -wjv . (19) 

Note that x T x = 1 implies wjw^- = 1 for each of the 8 update rules. The parameter < a < 1 controls the 
learning process. If j' and s' denote the values of j and Sj that minimizes the cost Eq. (19), the final step of the DLM 
algorithm is to set x = Wj> . 

In general, the behavior of the DLM defined by rules Eqs. (18) and (19) is difficult to analyze without the use of 
a computer. However, for fixed input messages yo = (2/0, 2/i) & n d yi = (2/2,2/3) for the and 1 event respectively, it 
is clear what the DLM will try to do: It will minimize the cost Eq. (19) by rotating its internal vector x to bring it 
as close as possible to y = (2/0,2/1,2/2,2/3)- After a number of events (depending on the initial value of x, the input y, 
and a), x will be close to y. However, the vector x does not converge to a limiting value because the DLM always 
changes its internal vector by a non-zero amount. It is not difficult to see (and supported by simulations, results not 
shown) that once x is close to y, it will keep oscillating about y = (2/0,2/1,2/2,2/3)- 

Let us denote by no the number of times the DLM selects update rule j = (see Eq.(18)). Writing Wq = 
(xq + S) 2 = 1 — a 2 + o?x\ and assuming that < q < 1, we find that the variable xq changes by an amount 
5 w (1 — a 2 )(l — Xq)/2x (neglecting terms of order S 2 ). If n is the total number of events then n — n is the number 
of times the DLM selects update rules j = 1,2,3. For j = 1,2,3 we have w 2 j = (x + S') 2 = o?x\ and hence x 
changes by 5' w —(1 — a 2 )x$/2. If x oscillates about y then xq also oscillates about 2/0- This implies that the number 
of times Xo increases times the increment must approximately be equal to the number of times x decreases times the 
decrement. In other words, we must have noS + (n — n )6' w 0. As x w y we conclude that n /n rts y 2 . Applying 
the same reasoning for the cases where the DLM selects update rule j — 1 shows that the number of times the DLM 
will apply update rules j = 0, 1 is proportional to j/g + 2/i- Therefore, the rate with which the DLM selects update 
rules j = 0,1 corresponds to the probability for observing a event in the quantum mechanical system. In other 
words, the DLM generates (1) events in a deterministic manner, with a rate that is proportional to the probability 
Po (Pi = 1 — Po) for observing a (1) event in the corresponding quantum mechanical system. Thus, a DLM is a 
simple "classical" dynamical system that exhibits behavior that is usually attributed to quantum systems. 

B. Stochastic learning machines 

The sequence of events that is generated by a DLM (network) is strictly deterministic. We now describe a simple 
modification that turns a DLM into a stochastic learning machine (SLM). The term stochastic does not refer to the 
learning process but to the method that is used to select the output channel that carries the outgoing message. 

In the stationary regime, the components of the internal vector represent the probability amplitudes. Comparing 
the (sums of) squares of these amplitudes with a uniform random number < r < 1 gives the probability for sending 
the message over the corresponding output channel. For instance, in the case of the Hadamard gate (see Fig. 5) we 
replace the DLM 2 by a SLM. This SLM generates a event if x\ + x\ < r and a 1 event otherwise. Although 
the learning process of this processor is still deterministic, in the stationary regime the output events are randomly 
distributed over the two possibilities. Of course, the rate at which different output events are generated is the same as 
that of the original DLM-network. Replacing DLMs by SLMs in a DLM-network changes the order in which messages 
are being processed by the network but leaves the content of the messages intact. Therefore, in the stationary regime, 
the distribution of messages over the outputs of the SLM-network is essentially the same as that of the original DLM 
network. 



IV. EVENT-BASED SIMULATION OF QUANTUM COMPUTERS 

A. Hadamard operation 34 

As an example of a DLM-based processor that performs singlc-qubit operations we consider the diagram shown in 
Fig. 5 (left). The presence of a message is indicated by an arrow on the corresponding line. The first component, called 
front-end, consists of a DLM that "learns" about the occurrence of and 1 events, meaning that the corresponding 
qubit is or 1, respectively. The second component transforms the data stored in the front-end and feeds this data 
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FIG. 5: Left: Diagram of the network of two DLMs that performs a deterministic simulation of a Hadamard gate on an event- 
by-event basis. The arrows on the solid lines represent the input and output events. Dashed lines indicate the flow of data 
within the DLM-based processor. Right: Simulation results for the Hadamard gate shown on the left. The input events are 
either of type with message (cos ipo, sini/>o) or of type 1 with message (cos ipi, sin tpi). A uniform random number is used to 
generate the type of input events. The probability for a (1) event is po (pi = 1 —po)- Each data point represents a simulation 
of 10000 events. After each set of 10000 events, a uniform random number in the range [0, 360] is used to choose the angles 
tpo and tpi . Markers give the simulation results for the normalized intensity in output channel as a function of <f> = ipo — ipi ■ 
Open circles: po = 1; Bullets: po — 0.5; Open squares: po = 0.25. Lines represent the results of quantum theory (see Eq.(23)). 



into a second DLM called back-end. The back-end "learns" this data. The learning process itself is used to determine 
whether the back-end responds to the input event by sending out either a or a 1 event. None of these components 
makes use of random numbers. 

An event corresponds to the arrival of a particle in either the or 1 state. The message is a unit vector y = (ya, y\) 
of two real numbers. We denote the number of (1) events by Nq (Aq) and the total number of events by A^ = Nq+N\. 
The correspondence with the quantum system is rather obvious: the probability for a event is given by |an| 2 w N /N 
and yo = Re ao/|ao| and y\ = Im a /|a |. The probability for a 1 event is Ni/N rts \a\\ 2 and y 2 — Re Oi/|oi| and 
y 3 = Im ai/|ai|. 

From Fig. 5 (left) it is clear that the transformation is just the real-valued version of the complex-valued matrix- 
vector operation that corresponds to the Hadamard gate (see Eq.(9)). A processor that performs the general singlc- 
qubit operation Eq.(8) is identical to the one shown in Fig. 5 (left) except for the transformation stage. For instance, 
to implement the X operation (see Eq.(9)) we only have to replace the transformation matrix of the Hadamard 
operation 




(20) 



by 








(21) 



In Fig. 5 (right) we present results of simulations using the processor depicted in Fig. 5 (left). Before the first 
simulation starts we use uniform random numbers to initialize the two four-dimensional internal vectors of DLM 1 
and DLM 2. Each data point in Fig. 5 (right) represents a simulation of 10000 events. All these simulations were 
carried out with a = 0.99. For each set of 10000 events, two uniform random numbers in the range [0, 360] determine 
two angles ipo an d V'l that are used as the message yo = (cosf/'O) sin^ ) (yi = (cos^i, sin ipi)) for the input event of 
type (1). Uniform random numbers are used to generate (1) input events with probability p$ (pi = 1 — po)- This 
corresponds to the input amplitudes ao = y/poe 1 ^ and a\ = pie 1 ^ 1 in the quantum mechanical system. 
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FIG. 6: Simulation results for the DLM-implementation of the network shown in Fig. 2. The input are events with message 
(cosi/>o,sini/>o). A uniform random number in the range [0,360] is used to choose the angle ipo- Each data point represents 
10000 events (No + TVi = N2 + N3 = 10000). Initially the rotation angle <f> — and after each set of 10000 events, <f> is increased 
by 10°. Markers give the simulation results for the normalized intensities as a function of (j>. Open squares: N /(N + Ni) (0 
events); Solid squares: N2KN2 + N3) (0 events); Bullets: Ns/(N2 + N3) (1 events). Lines represent the results of quantum 
theory. 



According to quantum theory, the probability amplitude b (b\) for the (1) output event is given by 



_Woo + ai\ = J_/ 1 1 
y/2 \ 01 -ao J V2 V 1 ~ X 



a 
«i 



(22) 



As the DLM-based Hadamard gate generates Nq/N events of type and Ni/N events of type 1, it is obvious from 
Fig. 5 (right) that these ratios are in excellent agreement with the probabilities 



1 + 2^p (l -Po) cos(?/>o - Vi) 



\h\ z = 



1 - 2- v /po(i -Po) cos (-00 ~ ^i) 



(23) 



as obtained from Eq.(22). 



B. Mach-Zehnder interferometer 34 



As a second example of event-based simulation of single-qubit operations we consider the Mach-Zehnder interfer- 
ometer network shown in Fig. 2. We use the equivalent DLM-based processor for the X operation. The phase-shift 
operation R((f>) is carried out by a passive device (that is, a device without DLMs) that simply passes messages of 
events and transforms messages of 1 events by performing a plane rotation about cj) of the two-dimensional vector 
representing the message. Thus, R{4>) transforms the message y = (yo,2/i) carried by a 1 event according to 



2/0 



1 
71 



COS( 



- smc 

COS( 



Uo 
Vi 



(24) 



In Fig. 6 we present a few typical simulation results for the Mach-Zehnder interferometer built from DLMs. We 
assume that the input receives events only, each event carrying the message (cos tpo> sin^o)- This corresponds to 
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FIG. 7: Diagram of a DLM-based processor that simulates a CNOT gate on an event-by-event basis. 



(ao,a\) = (cos -00 + isimpo,0) in the quantum system. We use a uniform random number to determine ipo. In all 
these simulations a = 0.99. The number of (1) events after the first X operation is denoted by No {N\). The 
number of (1) events after the last X operation is denoted by 7V 2 (N 3 ). The data points in Fig. 6 are the simulation 
results for the normalized intensity N/ (No + N\) for i = 0, 2, 3 as a function of fa Lines represent the corresponding 
results of quantum theory (see Eq.(12)). From Fig. 6 it is clear that the event-based processing by the DLM network 
reproduces the probability distribution as obtained from Eq.(12). Messages generated by DLMs preserve the phase 
information that is essential for the system to exhibit quantum interference effects. 

Summarizing: We have shown that a DLM-network can simulate single-photon quantum interference particle-by- 
particle without using quantum theory. In practicular, the previous example demonstrates that locally-connected 
networks of processing units with a primitive learning capability are sufficient to simulate, event-by-event, the single- 
photon beam splitter and Mach-Zehnder interferometer experiments of Grangier et al. 13 . The parts of the processing 
units and network map one-to-one on the physical parts of the experimental setup and only simple geometry is used 
to construct the simulation algorithm. In this sense, the simulation approach we propose satisfies Einstein's criteria 
of realism and causality 8 . 



C. CNOT operation 

The schematic diagram of the DLM-network that performs the CNOT operation on an event-by-event (particle-by- 
particle) basis is shown in Fig. 7. Conceptually the structure of this network is the same as in the case of a system of 
a single qubit. As input to the DLM-network we now have four (0,1,2 or 3) instead of two different types of events, 
corresponding to the quantum states |00), |01), |10), |11). Each event carries a message consisting of two real numbers 
y = (cos fa, sin fa) for i = 0, . . . , 3, corresponding to the phase of the quantum mechanical probability amplitudes, 
that is ao /| ao |, ■ • ■ , a 3 /|a 3 |. The internal state of each DLM is represented by a unit vector of eight real numbers 
x = (x , . . . , xr) and there are sixteen candidate update rules ({j = 0, . . . 7; Sj = ±1}, see Eq.(18)) to choose from. 
The rule that is actually used is determined by minimizing the cost function C = — sjwjv . The transformation 
stage is extremely simple: According to Eq.(13), all it has to do is swap the two pairs of elements (£2,23) and (xq,Xt). 

Instead of presenting results that show that the DLM-processor of Fig. 7 correctly performs the CNOT operation on 
an event-by-event basis, we consider the more complicated network of four Hadamard gates and one CNOT gate shown 
in Fig. 3. Quantum mechanically, this network acts as a CNOT gate in which the role of control- and target qubit have 
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TABLE II: Simulation results for the DLM-network shown in Fig. 3, demonstrating that the network reproduces the results 
of the corresponding quantum circuit, i.e. a CNOT operation in which qubit 2 is the control qubit and qubit 1 is the target 
qubit 6 . The first half of the events are discarded in the calculation of the frequency ft for observing an output event of type 
i = 0, 1, 2, 3, corresponding to the probability to observe the quantum system in the state 1 00) , |01), |10), or |11), respectively. 
For 200 events or more, the difference between the event-based simulation results and the corresponding quantum mechanical 
probabilities is less than 1%. 
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TABLE III: Simulation results for the DLM-network shown in 


Fig. 3 in which each back-end DLM of the individual gates has 


been replaced by a 


SLM. The latter uses random numbers to randomize the order 


in which different output events are generated 


but does not chani 


je the frequencies of the events. The first half of the input events are discarded for the calculation of the 


frequencies /; for observing an output event of type i = 0, 1, 


2, 3, corresponding to the probability to observe the quantum 


system in the state 00), |01), 10), or 


|11), respectively. 
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1 


0.010 


0.008 
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0.966 


2000 
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1 


0.005 


0.016 


0.963 


0.016 



been interchanged 6 . For this DLM-network to perform correctly it is essential that the event-based simulation mimics 
the quantum interference (generated by the Hadamard gates) correctly. In Table II we present simulation results for 
the DLM-network shown in Fig. 3. Before the simulation starts, we use uniform random numbers to initialize the 
internal vectors of the DLMs (ten vectors in total). For simplicity, we take as input messages 4>q = §\ = (f>2 = <fe = 0. 
All these simulations were carried out with a = 0.99. From Table II it is clear that, also for a modest number of input 
events, the network reproduces the results of the corresponding quantum circuit, i.e. a CNOT operation in which 
qubit 2 is the control qubit and qubit 1 is the target qubit 6 . 

As an illustration of the use of SLMs, we replace all the back-end DLMs in the CNOT circuit shown in Fig. 7 by 
SLMs. and repeat the simulations that yield the data in Table II. From Table III we conclude that the randomized 
version generates the correct results but significantly more events are needed to achieve similar accuracy as in the 
fully deterministic simulation. 

D. Number factoring 

Finally, we discuss the results obtained by a DLM-based simulation of the number factoring circuits depicted in 
Fig. 4. DLM networks (not shown) that perform the Toffoli gate operation and the Fourier transform, which involves 
several Hadamard operations and controlled phase shifts, are readily constructed by mimicking the procedure for the 
construction of the DLM network of the CNOT gate. Having done this, building the circuit in Fig. 4 entails nothing 
than connecting the DLM networks that simulate the various quantum gates. As this circuit involves seven qubits, 
the internal vectors of the DLMs have 256 elements. 

Section HE shows that if a = 7, we expect to find Qi = (Qf) — 0, Q 2 — (QZ) = 0.5, and Q 3 = (Qf) = 0.5. 
Similarly, for a = 11 wc expect to find Qi = (Qf) = 0, Qi = (Qf) = 0, and Qz = (Q§) = 0.5. In the DLM approach, 
simply counting the number of 1 events in the three output channels of the Fourier transform (inside the dashed box 
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FIG. 8: Event-by-event simulation of Shor's quantum algorithm for factoring the integer N = 15, using the value a = 7 (see 
Section HE). Each data point represents the average of 100 output events. The parameter that controls the learning process 
of the DLMs is a = 0.99. Left: Deterministic simulation employing DLMs. Right: Stochastic simulation employing SLMs. 



in Fig. 4) and dividing these numbers by the total number of events analyzed, we obtain numerical estimates for the 
qubits Qi = (Ql), Q2 = (Q2); an d Q3 = (QVl- ^ n Fig s - 8-11 we present simulation results for the DLM (left panel) 
and SLM (right panel) implementation of the circuit depicted in Fig. 4, for the two cases a = 7, 11 and for two choices 
of the control parameter a = 0.99, 0.999. After processing a few events, (less than 200 if a = 0.99 less than 2000 if 
a = 0.999) all the DLM networks reproduce the results of quantum theory with high accuracy. Replacing DLMs by 
their stochastic equivalents (SLMs), we see that the fluctuations are larger than if we use DLMs and that the system 
as a whole "learns" much slower. This is to be expected: the probabilistic mechanism to distribute events over the 
output channels of the SLMs makes much more "mistakes" than the deterministic process used by the DLMs. DLMs 
encode the information about the probability and phase in a much more effective, compact manner than SLMs. In 
the case of the latter, the correct probability distribution is encoded in a statistical manner and can only be recovered 
by analyzing a lot of events. 

From the description of the learning process, it is clear that a controls the rate of learning or, equivalently, the 
rate at which learned information can be forgotten. Furthermore it is evident that the difference between a constant 
input to a DLM and the learned value of its internal variable cannot be smaller than 1 — a. In other words, a also 
limits the precision with which the internal variable can represent a sequence of constant input values. On the other 
hand, the number of events has to balance the rate at which the DLM can forget a learned input value. The smaller 
1 — a is, the larger the number of events has to be for the DLM to adapt to changes in the input data. The results 
depicted in Figs. 8-9 and Figs. 10-11 confirm this behavior. 



V. DISCUSSION 



We have shown that locally connected networks of machines that have primitive learning capabilities can be used 
to perform a deterministic, event-based simulation of quantum computation. On the other hand it is known that the 
time evolution of the wave function of a quantum system can be simulated on a quantum computer 6,35 . Therefore, it 
is possible to simulate real-time quantum dynamics through a deterministic event-based simulation by constructing 
appropriate DLM-networks. The work presented in this paper suggests that there exist deterministic, particle-like 
processes that reproduce quantum mechanical behavior. 

Just as any other method for simulating quantum computers 30 , the DLM-based simulation approach requires 
memory resources that increase exponentially with the number of qubits. This exponential increase is merely a 
combinatorial effect and is not at all related to the quantum nature of the phenomena we want to simulate. As 
a matter of fact, it is present in all classical or quantum many-body systems (including quantum computers). For 
instance, if we consider one of the most basic statistical mechanics models, the Ising model, the number of possible 
states of this system also grows exponentially with the number of spins but there obviously is nothing "quantum" 
about this model. 
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FIG. 9: Same as Fig. 8, except that a = 11 (see Section HE). 
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FIG. 10: Event-by-event simulation of Shor's quantum algorithm for factoring the integer N = 15, using the value a — 7 (see 
Section HE). Each data point represents the average of 100 output events. The parameter that controls the learning process 
of the DLMs is a — 0.999. Left: Deterministic simulation employing DLMs. Right: Stochastic simulation employing SLMs. 



The computational efficiency of the event-based approach is lower than the efficiency of algorithms that directly 
compute the product of the unitary matrices. This is hardly a surprise: The former approach simulates quantum 
behavior by generating individual events. The latter can only simulate the outcome of (infinitely) many of such events 
and provides no information about individual events 8 . An analogy may be helpful to understand the conceptual 
difference between these two approaches. It is well known that an ensemble of simple, symmetric random walks may 
be approximated by a diffusion equation (for vanishing lattice spacing and time step). Also here we have two options. 
If we are interested in individual events, we have no other choice than to simulate the discrete random walk. However, 
if we want to study the behavior of many random walkers, it is computationally much more efficient to solve the 
corresponding diffusion problem. 

Our event-based approach can be extended to mimic the effects of decoherence. In quantum theory, decoherence 
causes the loss of phase coherence 36 . In the DLMs that we describe in Section IV, a single parameter (a) controls the 
loss of memory, of both the probability and the phase. A simple extension would be to control the learning process 
of the probability and phase separately, using two control parameters. We leave this topic for future research. 
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FIG. 11: Same as Fig. 10, except that a = 11 (see Section HE). 
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